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Abstract 

The  basic  direction  underlying  this  work  is  the  continuing 
study  of  certain  fundamental  features  associated  with  the  nonlinear 
wave  propagation  which  are  motivated  and  arise  in  physical  problems. 
The  usefulness  of  the  work  is  attested  to  by  the  varied  applications, 
and  wide  areas  of  interest  in  physics,  engineering  and  mathematics. 

The  work  accomplished  relates  to  fluid  mechanics,  nonlinear  optics, 
multidimensional  solitons,  Painlev^  equations,  long  time  asymptotic 
solutions,  new  linearizations  of  the  KdV  equation,  Riemann- 
Hilbert  problems  etc. 

il)  Research  Objectives  and  Accomplishments 

The  continuing  theme  of  the  work  performed  under 
this  grauit  has  been  the  study  of  nonlinear  wave  propagation  associated 
with  physically  significant  systems.  The  work  has  important  appli¬ 
cations  in  fluid  dynamics  (e.g.  long  waves  in  stratified  fluids) ,  non¬ 
linear  optics  (e.g.  self-induced  tranparency,  and  self-focussing  of 
light) ,  and  mathematical  physics  as  well  as  important  consequences 
in  mathematics.  Individuals  working  with  me  and  hence  partially 
associated  with  this  grant  include:  Dr.  i  :sios  Fokas,  Assistant 
Professor  of  Mathematics  and  Computer  Science.,  Dr.  Akira  Nedcaroura, 
post  doctoral  fellow  in  Mathematics  and  Computer  Science,  Mr.  Thied? 
T£Uia  and  Mr.  Paulo  Santini,  graduate  students  in  Mathematics  and 
Computer  Science.  Attached  please  find  the  technical  section  of 
our  recent  proposal  to  A.F.O.S.R.  In  this  proposal  many  of  the 
main  directions  and  results  are  outlined.  In  addition  Vv^e  nave 
listed  all  our  publications  and  preprints  carried  out  emd  conceived 
during  the  period  of  this  grant. 
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i'iniilly  \ie  point  out  an  entirely  new  and  important  area 
of  study  which  was  developed  during  the  past  year.  Namely  we 
have  employed  a  novel  linear  integral  equation  wiiich,  in  principal 
allows  one  to  capture  a  far  larger  class  of  solutions  to  the 
Korteweg-deVries  equation  (KdV)  than  does  the  v/ell  known  Gel'fand- 
Levitan-Marchenko  equation.  In  our  recent  paper  attached  we 
(a)  give  a  direct  proof  of  the  above  facts;  (b)  discuss  hov;  the 
Gel ' f and-Levitan-Marchenko  equation  may  be  recovered  as  a  special 
case,  (c)  characterize  a  three  parameter  family  of  solutions  to 
the  self-similar  O.D.E.  associated  with  KdV  and  which  may  be 
directly  related  to  the  second  Painleve  equation.  (We  note  that 
the  Gel ' fand-Levitan' Marchenko  equation  only  characterizes  a  one 
pareuneter  family  of  such  solutions) .  in  order  to  carry  out  (c)  we 
had  to  investigate  a  concrete  singular  integral  equation  in  which 
the  contour  L  consists  of  5  rays  all  passing  thru  the  origin.  The 
analysis  requires  the  full  power  (and  some  extensions)  of  the 
classical  theory  of  singular  integral  equations. 

It  should  be  remarked  that  (i)  the  integral  equation 
applies  to  potentials  of  the  Schrodinger  equation,  even  without 
the  application  to  KdV  or  (ii)  the  motivation  for  developing 

such  eui  integral  equation  originates  from  the  concept  of  summing 
perturbation  series.  (iii)  Recently  Flasclika  and  Newell  considered 
Pjj  via  monodromy  theory.  In  their  work  they  derive  a  formal 
system  of  linear  singular  integral  equations  for  the  general  solution 
of  Pjj*  However  the  highly  nontrivial  question  of  existence  of 
solutions  was  left  open.  How  their  work  and  ours  relate  is  a 
question  for  future  research.  (iv)  The  linear  version  of  KdV; 
u^+Uj^^j^  ®  0  is  solved  in  full  generality  as  a  special  case. 
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Some  future  directions  are;  (a)  Investigation  of  the  full 
generality  of  the  solutions  of  KdV  via  this  new  formulation. 

(b)  Development  of  sirailar  types  of  integral  equations  for  other 
nonlinear  evolution  equations,  as  well  as  ones  which  relate  to 
natural  "equilibrium"  states  for  KdV,  other  than  the  zero  (or 
vacuum)  state 

The  other  areas  of  interest  are  continuing.  They  include: 

.  Investigation  of  a  class  of  physically  significant 
nonlinear  singular  integro-dif ferential  equations,  and  associated 
novel  scattering  problems. 

.  Transverse  Instability  of  One-Dimens ioi.al  Transparent 
Optical  Pulses  in  Resonant  Media. 

.  Perturbations  of  Solitons  and  Solitary  Waves. 

.  Focussing  Phenomena  in  Nonlinear  Wave  Propagation. 

.  Two-dimensional  Lumps  and  Multidimensional  I.S.T. 

.  A  Connection  between  Nonlinear  Evolution  Equations 
and  Nonlinear  O.D.E.'s  of  Painlev^  Type. 


.  Discrete  I.S.T.  and  Numerical  Applications. 
.  Long  Time  Asymptotic  Solutions. 


.  Applications  of  Hirota's  Bilinear  Theory. 
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0 .  Forward 

The  main  purpose  for  the  continuation  of  this  research 
funding  is  to  support  the  work  presently  being  carried  out  by 
Mark  Ablowitz,  and  his  associates,  in  the  Mathematics  and  Computer 
Science  Department  at  Clarkson  College  of  Technology.  The 
principal  investigator  has  been  v;orking  in  the  general  area  of 
nonlinear  wave  propagation  for  almost  ten  years.  The  scope  of 
the  work  is  broad,  although  it  has  as  its  principal  focus  the 
uhderstudying  of  nonlinear  phenomena  connected  v;ith  the  v/ave 
propagation  which  arise  in  physical  problems.  In  recent  years 
significant  breakthroughs  have  been  made  and  this  area  of  research 
is  of  current  interest  to  mathematicians,  physicists,  and  engineers 
alike\  During  the  past  year,  the  active  research  funds  allowed 
us  to  support  Dr.  J.  Satsuma  as  a  Postdoctoral  Research  Associate 
at  Clarkson,  while  he  was  on  leave  from  Kyoto  University.  Dr. 

Satsuma  is  a  w’ell  knov:n  authority  in  this  field  of  research  and 
he  has  been  a  valuable  asset  to  our  research  programi.  Dr.  Satsuma 
returned  to  Kyoto  on  September  30,  1980;  nevertheless  we  have 
continued  to  collaborate  on  problems  of  joint  interest.  Subsequently 
Dr.  Akira  Nakamura,  also  from  Kyoto  has  come  to  Clarkson  as  a 
Postdoctoral  Associate  on  this  grant.  Dr.  Nakamura  has  written 
a  number  of  significant  papers  in  this  area  and  we  expect 
collaboration  between  us  to  be  very  fruitful. 

The  proposal  is  divided  as  follows.  In  the  first 
section  an  abstract  of  the  research  is  given.  In  the  second  section 
we  give  a  report  p.f  current  and  proposed  research.  The  third 
section  gives  references;  the  fourth  section  contains  curriculum 
vitaes  of  the  principal  investigator  and  Dr.  Nakamura,  and  the  fifth 
section  contains  a  proposed  budget  for  two  years. 


,\ 
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1 .  Abstrtict 

During  the  past  tv.’o  decades  significant  advances  in  t;;e 
study  of  nonlinear  v;ave  pher.or.ena  have  occurred.  These  advances 
have  allowed  researchers  to  begin  to  understand  some  of  th.c 
fundamental  building  blocks  associated  v;ith  nonlinear  waves  as 
well  as  being  able  to  obtain  solutions  to  a  number  of  nonlinear 
evolution  equations.  It  is  important  to  recognize  that  these 
studies  are  generic  in  nature  and  apply  to  numerous  physical 
problems  such  as  propagation  of  long  waves  in  stratified  fluids, 
self-f ocussing  in  nonlinear  optics,  self-induced  transparency, 
water  waves,  plasma  physics,  etc. 

In  the  period  of  time  mentioned  above,  both  approximate 
and  exact  methods  of  solution  to  problems  of  physical  sig¬ 
nificance  have  emerged.  Especially  significant  amongst  the 
exact  methods  of  analysis  is  v;hat  I  shall  refer  to  as  the 
Inverse  Scattering  Transform  and  the  associated  concept  of 
the  soliton.  This  method  has  found  applications  to  physics, 
engineering  and  mathematics  alike.  The  results  already  obtained, 
and  the  v;ids  ranging  interest  in  these  problem, s  have  motivated 
our  work.  In  this  proposal  we  discuss  some  of  the  research 
problems  v;hich  v;e  are  actively  pursuing. 
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2 .  Ciirrcnt  and  Proposed  P.osoarch 

Since  this  reserch  began  to  be  supported  by  the  Air  Force 
Office  of  Scientific  Research,  v;e  have  actively  studied  a  nur.'.ber 
of  problems  in  nonlinear  v;avc  theory.  In  v;hat  follov.’s  we  shall 
list  some  of  the  areas  which  we  have  studied  along  with  the 
principal  results  and  future  directions. 

(a)  A  Class  of  Physically  Significant  Singular  Nonlinear 
Integro-Dif f erential  Equations. 

Very  recently  we  have  become  interested  in  a  class  of 
nonlinear  singular  integro-diff erential  equations.  The  particular 
physical  application  of  this  problem  is  to  long  internal  gravity 
waves  in  a  stratified  fluid.  Hov;ever  both  the  v/ay  in  which  it 
arises,  and  the  relevant  mathematics  strongly  suggest  that  many 
other  applications  will  be  found  as  well.  In  fact  private 
communications  have  indicated  that  there  are  applications  to 
plasma  physics.  The  specific  equation  we  have  considered  is: 


where  T(u)  =  /  (-■ ^)  coth  u  (5)  d^. 


^  ^  represents  the  principal  value  integral  and  o  is  a  parameter , 
References  [1/2]  discuss  the  derivation  of  (2)  in  the  context 
of  internal  waves  •  As  6-^0  we  have  the  KdV  equation 


=  0, 


(2) 


wheras  if 


If*  < 


h  e  the  so-called  Bfenjamin-Ono  equation 
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u^+2uu^.+l!  -  0,  (3) 

v;here  H(u)  =  i/  is  the  Hilbert  transforn  of  u. 

Thus  equation  (1)  contains  as  limiting  forms  both  the 
KdV  and  Benjamin-Ono  equations.  The  fact  that  (1)  has  multi- 
soliton  solutions  ([3],  [4])  suggested  to  us  that  indeed  (1) 
may  be  solvable  by  the  Inverse  Scattering  Transform  (I.S.T.). 

In  our  recent  work  C5j,  [6]  v/e  have  found  a  Backlund  Transfor¬ 
mation,  a  generalized  Miura  Transformation,  soliton  and  rational 
solutions,  interesting  dynamical  systems  and  a  new  type  of 
scattering  problem.  This  scattering  problem  is  given  by  the 
equation 

iti/^+ (u- A)  =  u^~  (4) 

A 

+ 

where  u  satisfies  equation  (1)  ,  and  \li~  are  the  boundary  values 
of  a  function  analytic  in  the  strips  o<lmx<26  for  ij;'*’, 

-26<Imx<o  for  ip  ,  and  periodically  extended.  Specifically, 
equation  (4)  is  a  differential  Riemann-Hilbert  problem.  V.’hen 
A,y  are  given  by 

<lA  =  -kcoth2k6,  u  =  kcosech2k6, 

and  ip  (x)  =  ip^(x+2i6)  (by  periodicity)  we  find  that  in  the  limit 
6-*'0  we  have  the  Schrodinger  scattering  problem 


(5) 


v.’hich  is  the  linear  scattc-rin;;  probior'.  associated  v;iih 


•'uV  equat 


(2)  .  At  the  present  tine  v.-e  are  studying  the  inverse  scattering 
associated  v/ith  (4).  In  our  most  recent  v;ork  [7j  we  have  found 


Fredholm  integral  equations  for  the  Jost  functions  associated 
with  (4) .  This  is  unlike  previous  local  scattering  problems 
where  the  Jost  functions  satisfy  Voltcrra  type  integral  equations. 
V'e  note  that  nov;  the  scattering  problem.  (4)  and  associated 
nonlinear  evolution  equation  (1)  are  not  local.  In  those  cases 


where  the  Fredholm  equations  have  no  nontrivial  homogeneous 
solutions,  we  have  been  able  to  do  the  inverse  scattering  and  hence 
the  corresponding  initial  value  nrobJern  associated  v.’ith  (1).  This 
requires  both  6  and  max|u(x)t=0) 1  to  be  small  enough  (i.e., 
satisfy  certain  inequalities  -  in  some  sense  this  is  near  the 
KdV  equation)  .  When  6-*-oo  (the  Benjamin-Ono  limit)  we  have  found 
homogeneous  solutions  to  the  Fredholm  integral  equation.  Vve 
have  not  yet  carried  out  the  complete  inverse  scattering  analysis 
when  such  homogeneous  solutions  exist.  This  will  be  one 


important  aspect  of  our  future  work. 

Finally,  it  should  be  pointed  out  that  we  feel  that  there 
are  other  significant  nonlinear  singular  integro  differential 
evolution  equations  which  should  fall  into  the  category  of 
solvable  by  I.S.T.  and  processing  solitons,  we  shall  also 
investigate  such  possibilities  in  the  near  future. 

(b)  Transverse  Instability  of  One  Dimensional  Transparent 
Optical  Pulses  in  Resonant  Media. 

It  is  weir'known  that  ultrashort  optical  pulses  may 


propagate  coherently,  without  attenuation  in  certain  resonant 


media  CB,9].  This  phenomena  is  con-monly  rei erred  to  as  Se2:'- 
Induced  Transparency  (S.I.T.)  ana  has  been  intensively  studied 
experimentally,  numerically,  and  analytically  by  numerous 
researchers,  motivated  at  least  in  part,  by  significant  potential 
applications.  From  a  mathematical  point  of  view  the  one 
dimensional  equations  of  S.I.T.  are  very  special.  Namely,  it 
has  been  shown  that  these  equations  can  be  fully  integrated  by 
the  use  of  the  Inverse  Scattering  Transform  [10,111.  Specifically, 
the  above  analysis  has  shown  that  arbitrary  initial  values  break 
up  into  a  sequence  of  coherent  pulses,  which  do  not  decay  as 
they  propagate,  plus  radiation  which  rapidly  attentuates.  These 
coherent  pulses  are  referred  to  as  solitons. 

There  are  various  types  of  solitons  [8,91;  e.g.  "2ti 
pulses"  ("hyperbolic  secant  pulses"),  "on  pulses"  ("breathers") 
etc.  In  our  paper,  "Transverse  Instability  of  One-Dimensional 
Transparent  Optical  Pulses  in  Resonant  Media",  [12]  we  have 
shown  analytically,  that  the  2n  pulse  is,  in  fact,  unstable 
to  certain  transverse  variations  (i.e.  multidimensional  pertur¬ 
bations)  .  These  results  are  consistent  with  numerical  and 
experimental  studies  on  the  transverse  effects  in  S.I.T.  [13-141. 

The  latter  work  has  shown  that  transverse  variations  car  lead  -o 
frequency-amplitude  modulations,  and  in  some  cases  self-fccussing 
filaments.  Similarly  in  [151  we  have  recently  been  able  to  shov: 
that  the  breather  solution  (o~  pulse)  is  also  unstable  to  long 
transverse  perturbations.  Mathematically  speaking,  this  v.’ork 
was  difficult  because  the  earlier  analysis  had  to  be  much  further 
developed.  We  point  out  that  this  analytical  stability  calculation 
is  on  a  mode  which  is  much  more  complicated  than  a  permanent  ttavellin 
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viave  (i.e.  a  simple  soliton  2t!  pulse).  In  the  future  we  v;ish 
to  examine  the  stability  of  a  double  pole  solution  (i.e.  a 
limiting  form  of  a  breather  solution  just  before  it  breaks  up 
into  a  tv.’o  soliton  state)  as  well  as  attempting  to  more  fully 
understand  both  the  properties  of  the  tv;o  dim.ensional  model, 
as  well  as  looking  for  nev;  m.ultidimension  soliton  like  solitor.s 
(see  also  (e) ) . 

(c)  Perturbations  of  Solitons  and  Solitary  v.’aves . 

The  above  work  on  transverse  stability  of  solitons  in 
S.I.T.  led  us  naturally  to  the  problem  of  adding  general  weak 
perturbations  to  equations  which  admit  solitons  or  solitary 
waves  as  special  solutions  (both  in  one  and  more  than  one  dimension) . 
Some  of  the  mathematical  machinery  was  already  in  place  due  to 
the  work  done  in  part  (b)  described  above.  V’e  have  found  [16] 
that,  generally  speaking,  such  perturbation  problems  can  be 
successfully  handled  by  more  or  less  well  known  perturbation  methods. 
We  have  ccm.pared  our  results  to  some  of  those  in  the  literature 
which  employ  the  Inverse  Scattering  Transform  (see  for  example 
[17-19].  One  advantage  to  our  technique  is  that  it  also  applies  to 

problems  for  which  I.S.T.  does  not  apply. 

Our  analysis  shows  in  some  detail  that  there  is  quite 

different  phenomena  occurring  in  different  regions  of  space. 

Namely  near  the  peak  of  the  soliton  we  have  adiabatic  motion 
of  the  soliton  (or  solitary  wave) .  Away  from  the  soliton  a 
linear  W.K.B.  theory  applies.  The  results  are  asymptotically 
matched  in  order  to  obtain  a  uniformly  valid  theory.  To  our 
knowledge  this  theory  is  the  first  such  uniformly  valid  calculation 
of  a  perturbation  of  a  soliton  or  solitary  wave.  Previous  theories 
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,were  valid  in  linut:cU  ro:jic;‘.s  of  space  only. 

By  examinincj  other  equations  admitting  solitary  wave 
solutions  (i.e.  ones  v;hich  are  not  solvable  by  I.S.T.)  we 
believe  that  we  have  discovered  a  new  class  of  equations  v.’hich 
have  focussing  sngularities ,  (na-mely,  equations  which  have 
certain  solutions  which  are  "nice"  initially,  but  blow  up 
in  a  finite  time) . 

For  example,  we  have  discovered  evidence  that  strongly 
indicates  that  the  follov/ing  equation  is  in  this  class: 


u^+u^u  +u 

t  X 


XXX 


0 


(6) 


for  p V  4 .  V’le  hope  to  continue  to  investigate  such  questions 
in  the  near  future.  These  questions  are  of  both  mathematical 
and  physical  interest.  For  example,  such  a  question  arises  in 
the  propagation  of  water  waves  ( (d)  below) . 

(d)  Focussing  and  instability  associated  with  the 
propagation  of  water  waves. 

As  mentioned,  one  such  physical  problem  where  focussing 
occurs  is  that  of  water  waves.  We  have  considered  the  evolution 
of  gravity-capillary  v;aves  on  a  free  surface  of  a  layer  of  fluid 
with  constant  depth.  Beginning  with  the  standard  equations  of 
water  waves,  one  can  develop  the  evolution  equation  of  a 
modulated  weakly  nonlinear  periodic  wave  (with  fixed  central 
wave  number  k)  travelling  in  the  x-direction  whose  amplitude 
and  phase  vary  slov;ly  in  both  the  x  and  y  directions.  This 
problem  has  been  considered  by  Benney  and  Roskes  [20]  and  Davey 
and  Stewartson  [21]  without  surface  tension  and  by  Djordjevic 
and  Redekopp  [22]  who  included  this  effect.  For  sufficiently 
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?.on9  tin;es,  the*  complex  v.-uve  nr’plitucio  (x^y ,  t)  and  .moan  velocity 
potential  'J(x,y,t)  satisfy 


iA  +  o  ^  A  +  A 
t  1  :*:x  yy 


T  A 

X 


OL^ 


XX 


+  ^ 


yy 


=  0(1^1") 


(7) 


Here,  =  +1,  =  +l,a,g  are  real  constants  which  depend  on 

a  dimensionless  wavenumber  and  dimensionless  surface  tension 
coefficient.  Ablowitz  and  Segur  t23]  have  shov-’n  that  there 
exists  a  whole  parameter  regime  where  these  equations  have 
solutions  which  evolve  from  "nice"  initial  values  into  a  state 
where  the  amplitude  becomes  infinite  in  finite  time.'  In  this 
regime  the  fully  nonlinear  water  wave  equations  must  be 
analyzed.  Moreover,  in  [23]  it  was  demonstrated  that 
all  one  dimensional  soliton  solutions  to  (7)  are  unstable  with 
respect  to  long  transverse  variations.  These  results  indicate 
the  need  for  understanding  the  true  multidimensio.nal  character 
of  the  equation.  It  also  should  be  noted  that  the  v.'ork  in  [23] 
has  motivated  recent  experimental  work.  In  the  future  we  hope 
to  consider  hov;  these  focussing  solutions  can  be  understood  in 
the  context  of  the  fully  nonlinear  water  wave  equations,  as  well 
as  studying  potentially  solvable  multidimensional  cases  of  (7) 
(see  also  (e)  below) . 


(e)  Tv/o-dimcnsional  Lunps  in  Monlinear  Dispersive 


Systems 

Suitable  (long  v;ave)  limits  of  the  above  equation  of 
governing  watei:  waves  (7)  [21-22]  reduce  to 


o,?>  +  «>  =  -(  lA| 

1  XX  yy  '  '  X 


(8) 


02  =  +1/  or  for  very  long  waves  (comparable  to  the  v/eak 
nonlinearities)  : 


(u^  +  6uu  +  u  )  +  o^u  =  0/ 

'  t  X  XXX  X  3  yy 


(9) 


=  +1.  Both  of  these  equations  belong  to  tie  class  of  nonlinear 
evolution  equations  v/here  I.S.T.  is  applicable  [24,25J.  Indeed, 
recent  work  on  the  complete  integration  of  these  and  also  the 
three  wave  interaction  equations  has  been  undertaken  very  recently 
[see  for  example  26,27,281.  Alternatively,  via  Hirota's  method, 

N  plane  wave  soliton  solutions  can  be  directly  constructed  [see 
for  example  29] .  These  latter  solutions  are  quasi-one  dimensional 
they  describe  the  multiple  collision  of  N  solitons  each  of  which 
may  propagate  in  different  directions,  but  which  do  not  decay 
at  infinity. 

It  is  of  interest  to  find  essentially  two-dimensional 
solitons  which  would  be  analogues  of  those  in  one  dimension. 

By  taking  limits  of  the  N-plane  wave  soliton  solutions  described 
above  and  choosing  certain  parameters  appropriately,  we  have 
discovered  permanent  two-dimensional  nonsingular  lump  type 
solutions  decaying  in  all  directions  to  equation  (9)  [30, 
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sc“0  also  311.  Con'  in  liiiat.  (iiroction  wo  have  invoctin'-itecl 

analogous  solution:;  ociatec’  v-lth  ccjuation  (G)  .  VJe  have  shov/n 
that  eq.  (8)  possesses  a  lump  solution  of  envelope  hole  typo. 
Moreover,  v.’e  have  conntrvictcd  solutions  describing  multiple 
collisions  of  lumps  to  both  eqs.  (8),  (9)  [32].  It  is  of 

interest  to  see  how  tliese  solutions  can  be  fit  into  the  I.S.T. 
picture  and  also  to  investigate  other  equations  which  may  give 
rise  to  such  multidiir. -nsional  lump  type  solutions  (e.g.  S.I.T., 
the  self  dual  Yang  Mills  equations  etc.).  Me  hope  to  investigate 
such  questions  in  the  future  as  v;ell  as  to  study  the  generality 
of  the  initial  value  solutions  of  such  equations  obtained  by  I.S.T. 

(f )  A  Connection  Between  Nonlinear  Evolution  Equations 
and  Certain  Nonlinear  O.D.E.'s  of  Painleve  type 

•  The  development  of  the  inverse  scattering  transform  (I.S.T.) 
has  shown  that  certain  nonlinear  evolution  equations  possess 
a  number  of  remarkable  properties,  including  the  existence  of 
solitons,  an  infinite  set  of  conservation  laws,  an  explicit  sot 
of  action  angle  variables,  etc.  We  have  noted  in  [33]  that 
there  is  a  connection  betvjeen  these  nonlinear  partial  differential 
equations  (PDE's)  solvable  by  I.S.T.  and  nonlinear  ordinary 
differential  equations  (ODE's)  without  movable  critical  points. 

(Some  definitions:  a  critical  point  is  a  branch  point  or  an 
essential  singularity  in  the  solution  of  the  ODE.  It  is  movable 
if  its  location  in  the  complex  plane  depends  on  the  constants  of 
integration  of  the  ODE.  A  family  of  solutions  of  the  ODE  v/ithout 
movable  critical  points  has  the  P-property ;  here  P  stands  for 
Painlev6.)  In  [34-36]  we  have  announced  and  developed  a  number 
of  results  v/hich -indicate  that  this  connection  to  ODE's  of  P-type 
is  yet  another  remarkable  property  of  these  special  nonlinear  PDE's. 


/ 


-17- 

We  have  con joctured  thet: 

Every  nonlinear  ODE  obtained  by  a  similarity  reducticr; 
of  a  nonlinear  PDE  of  I.S.T.  class  is,  perhaps  after 
a  transformation  of  variables,  of  P-type. 

Here  we  refer  to  a  nonlinear  PDE  as  being  in  the  I.S.T.  class 

if  nontrivial  solutions  of  the  PDE  can  be  found  by  solving 

a  linear  integral  equation  of  the  Gel-fand-Levitar^-Marchenho 

form.  No  general  proof  of  this  conjecture  is  available  yet, 

but  we  have  proven  a  more  restricted  result  in  this  direction. 

It  is  known  that  under  scaling  transformations  certain  nonlinear 

PDE's  of  I.S.T.  class  reduce  to  ODE's.  Moreover,  the  solutions 

of  these  ODE's  may  be  obtained  by  solving  linear  integral  equations. 

We  have  shown  that  every  such  family  of  solutions  has  the 

P-property . 

Vfle  note  that  the  conjecture  in  its  strongest  form  relates 
to  ODE's  obtained  from  equations  solved  directly  by  I.S.T. 

There  are  many  examples  of  equations  solved  only  indirectly 
by  I.S.T.;  the  sine-Gordon  equation  is  one  of  the  best  knov:n 
examples.  An  ODE  obtained  from  an  equation  solved  indirectly 
by  I.S.T.  need  not  be  of  P-type,  but  it  may  be  related  through 
a  transformation  to  an  ODE  that  is. 

One  consequence  of  this  conjecture  is  an  explicit  test 
of  whether  or  not  a  given  PDE  may  be  of  I.S.T.  class;  namely, 
reduce  it  to  an  ODE,  and  determine  whether  the  ODE  is  of  P-type. 

To  this  end,  we  identify  certain  necessary  conditions  that 
an  ODE  must  satisfy  to  be  of  P-type  and  describe  an  explicit 
algorithm  to  determine  whether  an  ODE  meets  these  necessary 
conditions. 


Finally,  we  have  exploited  this  connection  in  order  to 
develop  both  solutions  and  asymptotic  connection  formulae  to 


some  of  the  classical  transcendents  of  Painleve  [  37]  as  v;cll 
as  others. 

In  the  future  we  shall  consider  the  follov;ing  problems: 

(1)  the  complete  connection  formulae  (i.e.  the  global  connection 
of  asymptotic  states)  for  the  interesting  Painleve  equations 
associated  v;ith  linear  Gel '  f  and-Levitan-Marchenko  equations; 

(2)  To  prove  that  the  O.D.E.'s  v;hich  v/e  have  derived,  in  fact 
satisfy  the  property  that  they  have  no  movable  essential 
singularities,  regardless  of  initial  conditions;  (3)  To 
develop  solutions  to  these  O.D.E.'s  which  correspond  to  general 
initial  conditions.  In  this  regard  the  recent  work  of  Flaschka 
and  Newell  [38j  may  be  of  interest.  (4)  Study  the  connection 
between  the  Backlund  transformations  developed  in  the  Russian 
literature  (see  the  review  [393)  and  by  Fokas  [403  and  their 
connection  to  l.S.T.  and  monodromy  preserving  deformations. 

(g)  Discrete  l.S.T.  and  Numerical  Schemes 

It  is  significant  that  many  of  the  concepts  related  to 
the  inverse  scattering  theory  apply  to  suitably  discretized 

nonlinear  evolution  equations;  for  example  the  Toda  lattice, 
and  discrete  nonlinear  Schrodinger  equation  (see  for  example 
[413,  [423).  It  is  of  interest  to  ask  whether  one  can  solve 
partial  difference  equations  (i.e.  numerical  schemes)  by  inverse 
scattering.  An  obvious  application  would  be  to  numerical 
simulations.  We  have  succeeded  in  analytically  devei-Oping  such 
schemes  [433.  These  schemes  can  be  shov;n  to  converge  to  a 
given  nonlinear  P.D.E.  (which  itself  is  solvable  by  inverse 
scattering)  in  the  continuous  limit.  Moreover  they  have  the 
nice  property  that  they  are  neutrally  stable,  have  exact  soliton 
solutions  and  possess  an  infinite  number  of  conserved  quantities. 
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Rocontly  v;e  have  (a  student,  Taha  and  myself)  begun  to 
compare  the  practical  numerical  simulation  of  a  given  nonlinear 
P.D.E.  (e.g.  cubic  nonlinear  Schrodinger  or  KdV)  using  traditional 
methods,  v;ith  our  newly  developed  schem.es.  In  this  sense,  v;o 
hope  to  assess  the  usefulness  of  various  numerical  schemes 
on  important  model  nonlinear  problems.  Preliminary  results 
indicate  very  encouraging  possibilities  for  these  "new” 
schemes  developed  via  I.S.T. 

(h)  Asymptotic  Solutions 

Aspects  of  the  asymptotic  solution  of  equations  solvable 
by  I.S.T.  have  been  discussed  by  many  authors.  For  a  review 
of  much  of  this  v;ork  the  reader  may  v\7ish  to  see  [37].  Despite 
all  of  the  work  already  done  on  this  question,  the  problem  of 
long  time  asymptotic  state  evolving  from  initial  data  containing 
both  solitons  and  the  dispersive  wavetrain  remained  unresolv'ed. 
Indeed  the  separate  questions  of  finding  the  asymptotic  states 
evolving  from  initial  conditions  containing  only  pure  solitons 
or  pure  dispersive  waves  had  been  solved.  The  difficulty  that 
must  be  overcome  when  both  states  are  present  is  that  the 
solitons  and  dispersive  wavetrain  are  of  differing  exponentially 
small  asymptotic  orders  in  certain  regions  of  space.  Recently 
[44]  we  were  able  to  completely  resolve  this  question  for  KdV. 
Moreover  we  find  explicit  formulae  giving  the  phase  shift  of  a 
soliton  v.’hen  it  interacts  ivith  both  solitons  as  well  as  the 
dispersive  waves.  A  corollary  to  this  result  is  the-def inition 
of  a  "perfect  soliton'  of  an  evolution  equation;  i.e.  one  which 
in  the  long  time  limit  interacts  elastically  with  any  sufficiently 
localized  disturbance. 


-20- 


V.’e  fool  Lht'  idc:ar-',  and  inetnoda  described  in  [44.1  r.lr;u2d 
apply  in  principle,  to  other  nonlinear  evolution  equations 
solvable  by  I.S.T.  V?o  hope  to  investigate  this  question  for 
some  of  the  other  nonlinear  evolution  equations  solvable  by 
I.S.T.  (e.g.  Nonlinear  Schrodinger,  Modified  KdV,  etc.) 

(i)  Recent  Work  by  A.  NAKA.MURA 

(a)  A  direct  method  of  obtaining  multiple  periodic 
wave  solutions 

Recently  the  study  of  multiply  periodic  wave  solutions 
to  nonlinear  evolution  equations  has  been  attracting  the  interest 
of  numerous  scientists.  Here  multiply  periodic  wave  solutions 
correspond  to  the  nonlinearly  superposed  state  of  several 
nonlinear  periodic  waves.  This  also  corresponds  to  the  periodic 
generalization  of  the  so-called  multiple  soliton  solutions  of 
nonlinear  evolution  equations. 

So  far  the  analytic  theory  to  calculate  these 
multiply  periodic  wave  solutions  has  been  based  on  the  theory  of 
rather  abstract  multidimensional  Riemann  surfaces.  Recently  v;e 
have  developed  a  completely  independent  and  different  approach 
to  this  problem,  based  on  Hirota's  direct  method  developed  in 
soliton  theory  [45].  Just  as  the  name  indicates,  our  method 
is  very  direct.  The  direct  method  is  constructed  by  elementary 
techniques.  By  this  direct  method,  v;e  have  obtained  multiply 
periodic  v?ave  solutions  [46]  of  the  intermediate  long  wave 
equation  [47]  which  includes  the  KdV  equation  and  Benjamin-Ono 
equation  as  limits.  This  intermediate  long  wave  equation  is 
considered  to  be- important  in  the  description  of  internal  wave 
soliton  pnenomena  and  has  recently  attracted  wide  interest  [48]. 
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At  prosouu,  our  direct  apriroccn  hov.’ovcr  has  boor;  developed 
to  calculate  only  to  a  2-poriodic  wave  solution.  V.’e  hope 
to  extend  our  direct  method  to  treat  arbitrary  h-pcriodic  v;ave 
solutions . 

(b)  Cylindrical  solitons,  solitons  in  nulti-diir.ensional 

systems 

Investigations  of  solitons  in  multi-dimensional  systems 
are  important  especially  in  view  of  the  actual  application  to 
the  real  physical  phenomena.  Here  just  as  the  KdV  equation 
has  crucial  importance  as  a  model  equation  in  one-dinensional 
nonlinear  dispersive  systems,  the  so-called  cylindrical  KdV  L49] 
and  spherical  KdV  [50]  equations  are  important  model  equations 
respectively  in  tv;o  and  three-dimensional  nonlinear  dispersive 
systems.  Both  cylindrical  and  spherical  solitons  have  been 
studied  experimentally  and  numerically  [49,51],  However  so 
far  exact  analytical  solutions  to  these  equations  are  either 
not  known  (in  the  case  of  spherical  KdV  equations)  or  could  not 
reproduce  experimental  and  numerical  results  (in  the  case  of 
cylindrical  KdV  equation  [53]).  By  generalizing  Hirota's  direct 
method  appropriately,  v;e  have  obtained  a  new  analytic  solution 
of  the  cylindrical  KdV  equation  which  does  correctly  reproduce 
experimental  and  numerical  results  [54].  We  have  also  calculated 
the  Backlund  transformation  for  the  cylindrical  KdV  equation  [55]. 
These  solitor.  solutions  have  a  certain  self  similar  character. 

Although  the  cylindrical  KdV  equation  has  been  successfully 
analyzed,  the  analogous  solutions  of  the  spherical  KdV  equation  is 
not  obtained  by  a  straightforward  application  of  the  method  to 
solve  the  cylindrical  KdV  equation.  However  since  the  cylindrical 
and  spherical  KdV  equations  are  very  similar  to  each  other  in  form, 
there  is  a  possibility  that  a  generalization  of  the  method  to 
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soivo  the  cylindrical  lUiV  eciuation  can  solve  the  spherical  one. 

We  are  considorinc;  such  a  generalization. 

c)  I'lultiple  decay  mode  solutions  of  nonlinear 
evolution  equations. 

In  the  case  of  the  cylindrical  KdV  equation  which 

has  been  discussed  above,  due  to  conservation  of  energy,  the 

inwardly  running  "ring-form"  soliton  grows  in  amplitude  and 

the  outv;ardly  running  one  decays.  Thus  it  is  natural  to  suspect 

that  the  similarity  type  solitons  which  appear  in  the  cylindrical 

KdV  equation  have  some  relationship  to  decaying  modes  of 

other  nonlinear  evolution  equations.  In  fact,  guided  by  this 

consideration,  very  recently  we  have  found  certain  similarity 

type  multiple  decay  mode  solutions  of  the  two-dimensional  KdV 
% 

(Kadomtsev-Petviashvili )  equation  [56],  Furthermore  it  has 
been  found  that  these  newly  found  decaying  mode  solutions  and 
usual  soliton  solutions  can  be  nonlinearly  superposed  with  each 
other.  Further  investigation  of  these  and  other  solutions  of 
two-dimensional  KdV,  and  their  relationship  to  I.S.T.  is 
under  consideration. 

d)  Chain  of  Biicklund  transformations 

It  is  well  known  that  the  typical  nonlinear  evolution 
equations  of  current  interest,  KdV  and  modified  KdV  (=mKdV) 
equations  are  related  v;ith  each  other  by  the  Miura  transform 
or  (non-auto)  Backlund  transform.  In  fact  this  process  continues 
to  other  equations  sequentially.  To  grasp  general  characteristics 
of  this  important  Backlund  transform,  it  is  worthy  to  study 
this  series  further.  We  have  recently  carried  out  the  study  of 
such  Backlund  transform  chains  e.g.  KdV-*-inKdV-»''' second  roKdV" 
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L  51  ]->■"  third  r.'.KdV"  .  .  .  [  58  J  whore  arrows  denote 
It  is  further  hoped  that  from  these  studies, 
to  extract  some  meaningful  information  about 
the  Backlund  transform  itself  L58J. 


j 

! 


Backlund  transforms, 
it  will  be  possible 
the  nature  of 
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form  way,  starting  from  an  arbitrary  classical 
Lie  group,  analytically  solvable  lattice  models, 
among  which  is  a  far*-reaching  generalization  of 
the  XXX  model,  given  by  a  quantum  Hamiltonian 
The  proposed  S-matrix  method  will  be 
instrumental  in  finding  energy  levels,  and  in  the 
computation  of  correlation  functions  for  the  known 
and  new  lattice  models  of  Toda  and XXX  type  con¬ 
sidered  here. 
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A  new  integral  equation  which  linearizes  the  Korteweg-de  Vries  and  Painleve  II  equa¬ 
tions,  and  is  related  to  the  potentials  of  the  Schrodinger  eigenvalue  problem,  is  present¬ 
ed.  This  equation  allows  one  to  capture  a  far  larger  class  of  solutions  than  the  GeTfand- 
Levitan  equation,  which  may  be  recovered  os  a  special  case.  As  an  application  this 
equation,  with  the  aid  of  the  classical  theory  of  singular  integral  equations,  yields  a 
three-parameter  family  of  solutions  to  the  self-similar  reduction  of  Korteweg  -de  Vries 
which  is  related  to  Painleve  11. 

PACS  numbers;  02.30. -t'g 

Since  the  work  of  Gardner  etaL  in  1967/  there  has  been  wide  interest  in  the  analysis  of  nonlinear 
evolution  equations  solvable  by  the  so-called  inverse-scattering  transform  (1ST).  The  prototype  exam¬ 
ple  is  the  Korteweg-de  Vries  (KdV)  equation 

(1) 

In  this  note  we  shall  present  a  new  linear  integral  equation  which,  in  principle,  allows  one  to  capture  a 
far  larger  class  of  solutions  than  does  the  Gerfand-Levitan  equation.  Specifically  we  claim  that  if 
^k;x,i)  solves 

V(k;  x,t)*i  exp[ i(kx*k»t )]  JT  dX(/ )  =  exp( , (kx (2) 

where  dX(Jb)  and  L  are  an  appropriate  measure  and  contour,  respectively,  then 

(3) 

solves  the  KdV  equation.  The  well-documented  physical  significance  of  the  KdV  equation,  of  its  self- 
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similar  analogue,  and  ol  the  associated  Scbrodingcr  scattering  problem  require  us  to  attempt  to  char¬ 
acterize  the  form  of  the  most  general  solution/potential  possible. 

We  now  enumerate  the  basic  results  given  in  this  note,  (i)  We  give  a  direct  proof  that  (2)  and  (3)  solve 
(1);  (ii)  we  show  how  the  well  known  Gel’fand-Levitan  equation  can  be  obtained  from  (2)  as  a  special 
case;  and  (iii)  we  characterize  by  a  matrix  Fredholm  equation  a  three-parameter  family  of  solutions  to 
the  similarity  ordinary  differential  equation  of  (1)  which  is  directly  related  to  the  classical  second  equa- 
tion'of  Painlevd  (P  n).  We  end  with  some  remarks  regarding  the  role  of  Backlund  transformations  and 
relevant  generalizations. 

We  now  consider  (i).  The  point  of  view  we  take  here  is,  in  spirit,  similar  to  that  of  Zakharov  and 
Shabat.*  Specifically,  by  direct  calculation  we  show  that  solutions  of  (2)  substituted  in  (3)  satisfy  (1). 

We  make  two  assumptions:  (a)  dX  and  L  are  such  that  differentiation  by  x,  t  may  be  interchanged  with 

(b)  the  homogeneous  integral  equation  has  only  the  zero  solution.  Defining  where  I, 

>  ♦  8,*,  after  some  manipulation  we  find 

L<((k;x,t)*ieKp[i(kx*k*t)]  f  ^^■dH,l)=3k[k(p^^-i<p^^^-iu<p].  (4) 

Similar  calculations  show  that  the  quantity  in  brackets  in  the  right-hand  side  of  (4)  satisfies  the  homo¬ 
geneous  integral  equation.  Hence  iufp^O  which  implies  £^^0,  whereupon  9gl2jiL<p)dX 

*0  is  (1).  Moreover  the  equation  Ai<p^0  is  directly  related  to  the  Schrodinger  eigenvalue  problem.  If 
we  define 

9{k;x,t)*^(k;x,  f )exp[f(*x+**<)/2j, 
then  gives 

(5) 

Next  we  pass  on  to  (ii).  The  classical  theory  of  inverse  scattering  and  appropriately  decaying  solutions 
of  KdV  may  be  most  easily  obtained  as  follows.  Let  the  measure  dX{k)  =  r,(ik)dk/2v,  where  r,(k)  is  the 
usual  reflection  coefficient  of  u(r,0)  and  the  contour  L  goes  over  all  the  poles  of  rj(,k).  [Here  we  have 
assumed,  for  convenience,  that  u(r  ,0)  —  0  rapidly  as  |  x|  —  ‘°.]  Then  substituting  the  expression  for  <p 
into  (2),  defining 

K{x,  y,  f  )=-(i)/x.^(*;  X,  f)exp[»(*y+A*t)/2]dX(*), 


and  using 

exp[  »■(*  +  ! )x/2]/{l ■*■  *)  a  -  f  { ji%xp(  Hb  + 1)^/2] /2di} 


(A,  I  satisfy  Imk,  Imf  >  0),  we  obtain 

K{x,y,t)*F[x*y;t)*  /f/f(x,  f ;  f  )F( «  +  y;  f  )d«  =  0, 
where 

F{x,  0 «(i)/r exp( i{kx/2*k»t )lrfA(*), 

and  u(x,  t)* 28,X(x, x;  t ).  Hence  by  choosing  the 
above  measure  dX  and  contour  L,  the  Gel'fand- 
Levltan  equation  (6)  may  now  be  completely  by¬ 
passed. 

Sollton  solutions  of  (1)  may  be  calculated  in  a 
particularly  easy  manner  from  (2).  Locations  of 
the  poles  on  the  imaginary  k  axis  in  r{k,  0)  corre¬ 
spond  to  sollton  amplitudes,  and  the  residues  of 
r{k,  0)  at  these  locations  play  the  role  of  the  nor¬ 
malization  coefficients.  Pure  solitons  may  also 
be  obtained  by  taking  the  measure  as 

dX(k),»ieMk-iK.)dk 

/■I 


(6) 


(L  passes  through  the  k*i Xy).  Then  (2)  reduces 
to  a  linear  algebraic  system  from  which  the  well 
known  N-soliton  solution  is  immediately  obtained. 

We  now  discuss  (iii).  The  KdV  equation  admits 
the  similarity  transformation  u(x, /)»0(x'y(30*^*, 
where  x'-x/(3l)'''*.  The  equation  for  U  is  given 
by  (dropping  the  primes) 

K,(£/)*i;''  +  6f/f/'-(2lf+x£f')»0.  (7) 

We  note  that  (7)  is  directly  related  to  PH: 

Pj(V)»  V*-xF-2V*»o.  (8) 

Specifically  we  note  that  the  transformations  U 
K*-  V,  F«(£;'+flt)/(2ff-x)relat8(8)tothe 
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equation 

20*  -  x£I+[i/ +  V' -  (in*]/(2t;  -  4  =>  0 

with  >'■0(0 -f  1).  However,  by  direct  calculation 
[(20-xVf,(0)l'«(20-x)i:,(0),  hence  K,(0)  is  an 
integral  of  (7),  and  thus  there  is  a  direct  trans¬ 
formation  between  (7)  and  (8).*  One  may  make 
use  of  these  transformations  fo  find  all  the  known 
(see,  for  example,  Lukashevich*  and  Erugln*) 
elementary  solutions  of  P  n.  Ablowitz  and  Scgur* 
had  established  a  connection  between  P II  and  1ST 
and  had  characterized  a  one-parameter  family  of 
solutions  via  the  Gel’fand-Levitan  equation.  Re¬ 
cently  Flaschka  and  Newell*  considered  P II  via 
monodromy  theory.  In  the  latter  work  the  authors 
derived  a  formal  system  of  linear  singular  inte¬ 
gral  equations  for  the  general  solution  of  P II. 
However,  the  highly  nontrivial  question  of  exis¬ 
tence  of  solutions  was  left  open. 

An  application  of  the  result  presented  above  in 
(i)  is  that  a  three-parameter  family  of  solutions 
of  (7)  may  be  obtained  from  the  linear  singular 
integral  equation 

tonL,  (9) 

where  d(0*/(^)“e*p{*(fjf+fV3)]  and 

(see  Fig.  1),  A>*Pa  =  Pi,  -^s  =  ^4  =  p2. 

(Hereafter y  always  stands  for  j  =  I, . . . , 
5).  The  solution  to  (7)  is  then  obtained  from 

{<p  depends  parametrically  on  x).  We  note  that 
both  (9)  and  U  are  obtained  from  (2)  and  (3)  by  a 


self -similar  reduction.  Moreover  the  contours 
Lj  are  obtained  by  finding  the  solution  to  the  lin¬ 
ear  problem  w''-(w  +  xw')*0  in  terms  of 

integral  representations  and  then  deforming  these 
contours  so  that  they  all  pass  throu^  the  origin. 
For  example,  note  that  L^*l^  may  be  deformed 
to  the  usual  Airy-function  contour.  If  we  restrict 
ourselves  to  this  Airy  contour,  tiie  result  in  Ref. 

5  is  obtained  in  the  same  manner  as  that  in  (ii) 
above. 

We  shall  proceed  to  demonstrate  that  (9)  may  be 
reduced  to  a  system  of  Riemann-Hilbert  problems 
which  are  sol^le  using  Fredholm  theory.  For 
this  we  need  the  full  power  of  the  classical  theory 
of  singular  integral  equations.^* 

Consider  the  sectionally  holomorphic  function 

“O’ 


(11) 

♦  *(!)=  toaLf, 

where  <^*(/)  for  t  on  Lf  has  the  standard  defini¬ 
tions^'*  of  limits  <rf  ♦(«)  as  *  -f  from  the  "left- 
hand  side’’  (+)  and  "right-hand  side"  (-)  of  Lj, 
and  where  principal-value  integrals  are  implied 
when  needed.  With  use  of  (11),  and  Eq.  (9)  for  t 
on  Lj  and  -/  on  -Lj,  we  obtain  a  system  which 
we  choose  to  write  in  the  form 

♦  *U)  =  C(/ )*-(/> +£(/),  tonZj,  (12) 

where  =  **(1)  ■[♦*(/),  ♦'(-!)]  *', 

r (0 = ♦  *(  -I ),  F(o  •  i7(/)/f(i ),  -/  ( -/  )H(  -1 )  ]  *• , 

if  t  on  Lj,  0  if  (  on  -Lj}  and  the  compo¬ 
nents  of  the  2x2  matrix  G(i)  are  G,,(l)  =  ~2b(f) 
xH(t)=-C„(-<),  Gu«C„«l. 

One  can  prove  the  following  statements. 

(a)  ♦*(-<),  ♦'(-<)  are  “minus"  and  "plus"  func¬ 
tions,  respectively,  (b)  Necessary  conditions  for 
solvability  of  (12)  are  the  symme^  conditions 
G(<)  «'lG(-<))"*,  £(f)+0(<)£(-f)"0,  which  are 
satisfied  by  the  above  G,  F.  (c)  Thus  (12)  defines 
a  system  of  discontinuous  Riemann-Hilbert  prob¬ 
lems  with  the  additional  restriction  that  4~(f ) 

However  this  conditien  can  be  relaxed 
since  one  can  show  that  (12)  always  admits  a 
solution  with  this  restriction,  and  moreover,  in 
our  case  the  solution  is  unique. 

In  order  to  solve  (12)  we  first  consider  the 


The  lines  of  discontinuity  of  4(z)  are  Lj-,  thus 
using  the  PlemelJ  formulas,  we  have 
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homogeneous  problem.  The  standard  procedure 
is  to  transform  the  discontinuous  homogeneous 
problem  to  a  continuous  one,  and  then  obtain  the 
fundamental  set  of  solutions. 

Associated  with  a  given  contour  define  the 
following  auxilliary  functions 


where  Zj  is  some  j  -dependent  fixed  point  off 
The  branches  of  the  above  functions  are  chosen 


such  that  ^  and  are  plus  and  minus  func¬ 
tions,  respectively  (c.g«,  the  branch  cut  for  ^ 
is  taken  between  0,  and  hence  lies  to  the 
right  of  Xj),  The  properties  (a^j|(0+)»expl-«rX^J, 
w^»(0-)»cxpliffX^J,  allow  us 

to  map  the  homogeneous  system  )♦“(/) 

which  has  a  discontinuity  at  /«0  to  the  following 
Riemann -Hilbert  system  which  is  continuous  at 
the  origin: 

/onl^,  (14) 

where  we  have  used  the  transformation  ) 

and  hence  ^/) 

=  lR^f)f‘A-‘C(7)An-(/),  with  A,  0*(f)  defined 
by 


exp[i»Xjj/2l  ^  expliirx^/2] 


expUsXji/2l 


e3q)U»X^/2l 


(15) 


where  for  J* 2, 3, 4  we  have  0*1,  ^*2  and  for  j  »1,  5  we  have  o  *2,  /3=sl;  the  X^ii  and  are  defined 

by 

exp[frX^i]=p^4(l+p^*)^*,  exp(ivXya]«-p^  ♦(l+p^*)*^*,  A^i  =  exp[2tiTX^j],  Aj,*exp[2tvX^,]« 


The  matrix ^(f)  has  the  properties  ^(f)  =  lg(-f)r*  and  detg  =  -l. 

One  may  characterize  a  solution  of  the  system  (14)  by  imposing  the  condition  *2^  as 

in  A*.  This  leads  to  a  Fredholm  equation  for,  say,  f^(/),  which  however  must  be  interpreted  In  a 
suitable  principal -value  sense  as  it  does  not  converge  in  the  normal  sense  at  infinity.  Alternatively, 
one  may  obtain  a  regular  Fredholm  equation  of  the  second  kind  by  imposing  conditions  at  a  finite  point 
off  all  contours,  say  a  » 1.  This  leads  to  the  following  Fredholm  equation  for 


t  on  Cj , 


(16) 


where  ^-*(1),  )»  and  ^  la 

the  unit  matrix.  Any  two  linearly  Independent  £ 
vectors,  say  di,a,  lead  to  a  fundamental  matrix 
j*(f)llor  the  system  (14). 

With  use  of  the  above  results  the  fundamental 
matrix  of  the  discontinuous  problem  (12)  is  given 

by 

X*(/)-An*(f)l4:, ♦(/),*,’(/)) .  (17) 

Hence  the  solution  of  (12)  is  given  by 

♦  *(/) 

*(/)/-  dr.  (18) 

Having  obtained  *\t)  and  using  (11)  to  obtain 
yK/),  we  have  characterized  a  three-parameter 
family  of  solutions  of  U.  With  use  of  the  results 
of  Fredholm’s  theory  the  nonmovable  critical- 


point  property  of  U  is  easily  verified. 

Finally,  we  make  some  remarks.  First,  we 
only  expect  from  (2)  to  obtain  solutions  to  P II  in 
the  range  -  j  <  a  <  3.  To  obtain  the  solution  for 
all  ranges  of  a,  we  believe,  the  Backlund  trans¬ 
formations  (following  Rosales‘S  and  *Tinite  per¬ 
turbations”  (see,  for  example,  Ablowitz  and 
Cornille*')  of  suitable  elementary  solutions  must 
be  employed.  Similarly,  wider  classes  of  solu¬ 
tions  to  KdV  should  be  obtainable  this  way  (we 
shall  remark  on  this  more  completely  in  the  fu¬ 
ture).  Second,  straightforward  generalizations 
to  the  higher-order  KdV  equations,  as  well  as  to 
many  other  nonlinear  evolution  equations,  are 
possible.  Third,  motivation  for  some  o(  the  ideas 
in  this  note  originate  from  the  concept  o£  sum¬ 
ming  perturbation  series.  Relevant  perturbation 
scries  can  be  readily  developed  (see,  for  exam¬ 
ple,  Refs.  lO^nd  12), 
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